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The progress for this period of the grant is best described in the 
attached papers and abstracts which- have been prepared. 

The first attachment is AIAA paper 85-0331, entitled "A Finite Volume 
Method for Calculation of Compressible Chemically Reacting Flows." It is 
to be presented at the AIAA 23rd Aerospace Sciences Meeting in Reno in 
January 1985. The previously reported algorithm for solving the stiff set 
of equations for chemically reacting flows is reported in detail. The 
computations have been do le for a two-dimensional geometry and, as we 
expected, show the. excellent computational efficiency obtained in the one- 
dimensional case. These calculations are for the Euler equations with the 
Rogers^Chinitrz hydrogen-oxygen model. 

The second attachment is an abstract submitted to the AIAA CFD Conference 
to be held in Cincinnati in July 1985. The finite volume method has been 
extended to viscous flows. Our results verify the accuracy of the non- 
reacting case by comparison with Carter's calculation. The reacting case 
for the attached boundary layer appears to be satisfactory, although there 
is no comparison case for validation. The calculations for recirculating 
flow are preliminary. 

The third~attachment is an abstract submitted to the AIAA 21st Joint 
Propulsion Conference, also in July 1985. It is the start of our study to 
investigate the 'effect of parameters such as fuel/air ratio, geometry, etc. 
on the physics of the flow. 

During the reporting period, Mr. Thomas Bussing spent about two months 
at NASA Langley Research Center during the summer of 1984. 



A FINITE VOLUME METHOD FOR THE NAVIER-STOKES 
EQUATIONS WITH CHEMICAL REACTION 

Thomas R. A. 3ussing* 
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ABSTRACT 


Over the past two decades Supersonic Combustion Ramjets (Scramjets) 
have received considerable attention; One critical element in the design 
of a. Scramjet is the de'tailed understanding of the complex flow field in 
the engine during various phases of operation. One area of interest is 
the computation of chemically. reacting flows in the vicinity of flame 
holders. The purpose of this proposed paper is to study the characteristics 
of a method for solving the Navier-Stokes equations with chemical reaction. 


The flow field in a Scramjet is governed by the Navier-Stokes equations 
coupled to the finite rate chemistry equations. These equations represent 
different processes which can have their own characteristic time scales. 

In the problems of interest the time scales can be orders of magnitude 
different (stiff) leading to severe time step restrictions on the compu- 
tational procedure. If only the steady state solution is desired several 
steps can be taken to accelerate the iterative solution. The details of 
these steps for the Euler equations with Finite Rate Chemistry are given 
in reference fl], The numerical procedure is based on a point implicit, 
version of ♦‘he Jameson, Schmidt and Turkel [2] finite volume/multistage 
scheme. 

Briefly the idea behind the point implicit method involves modifying 
the original stiff sot of equations to another set where the time stiff- 
ness associated with the chemical time scales has been removed. This can 
be accomplished by multiplying the unsteady terms by a scaling matrix. 

For example, instead of solving the stiff equation system given by 



we. could solve 


( 1 ) 


_ 3 jt -3x 

where 5 is picked to be the desired scaling matrix. A convenient choice 
for 5 is the matrix obtained by numerical lv treating the chemical source 
terms implicitly. Hence, the name point implicit. The technique trans- 
lates to solving the equations in pseudo time, marching each physical 
process at its own respective characteristic time scale. Note equations 
(1) and (2) have the same steady state solutions. 


( 2 ) 


The. method is .being applied to reacting boundary layer and rearward 
facing step flows. The full™ Navier-Stokes equations with H^-Air finite 
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rat': chemistry are used. 
Chmitz [3] and consists 


The chemistry mccel was dove I 
of the. following two steps: 


ted by Rogers and 


H 



2 OH 



+ 2 OH 


2 H 2 0 


Tentatively we expect to apply the method to chemically reacting flows 
over flat plates and rearward facing steps. A few preliminary examples of 
these types of flows will now be given.-. The first example compares the 
predicted non-reacting solution to a numerical data set generated by Carter 
[4J. The.. results are for supersonic flow over a flat plate. As the figure 
illustrates, our computed velocity profile agrees well with Carter's data. 

•If -chcnrcai re a c tore rPi'S" allowed to occur, then an H 0 species contour like 
that shown in Figure 2 is produced. Note as expected the reaction is 
largest next to the surface where the static temperature is highest. 

Finally a chemically reacting flow over a rearward facing step is shown in 
Figure 3 (M. =3, T. = 900K, T . ^ - 1QQ0K, -11 = 0.1). Each of these examples 

would be considered in greater 1 ^detai 1 in the proposed paper. With the 
method proposed in this abstract, convergence to steady state was achieved 
in approximately SCO iterations compared to approximately 5 x 10® for a 
purely explicit scheme. 


TABLE Or SYMBOLS . 

* . o 

HR J (Z Y. Kr.)dA/(Cv. T. + .5 tT) 
J i i i m - in 

$ S techier* trie mixture ratio 


o of srecios i/mixturo density 

l 


KF Heat of formation 
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ABSTRACT 


NUMERICAL INVESTIGATION OF 2-DIMENSICNAL H -AIR 
FLAME. HOLDING OVER RAMPS AND REARWARD FACING STEPS 

Thomas R. A. Bussing* 

Earll M. Murman* 

Massachusetts Institute of Technology 
Cambridge, Massachusetts 02139 

Over the past two decades Supersonic Combustion Ramjets (Scram jets) have 
received considerable attention. One critical element in the design of a 
Scram jet is the detailed understanding_of the complex flow field in the 
engine during various phases of operation. One area of interest is the com- 
putation of chemically reacting flows in the vicinity of flame holders. The 
purpose of this proposed paper is to study the flame holding characteristics 
of simple ramps and rearward facing steps. Both of these configurations are 
considered candidates for Scram jet flame holders. 

_ The flow field in_ a Scramjet is g over ned by the Navicr-Stokes equations 
coupled to the finite rate chemistry equations* These equations represent 
different processes which can have their own characteristic time scales. 

In the problems of interest the time scales can be orders of magnitude 
different (Stiff) leading to severe time step restrictions on the computational 
procedure. If only the steady state solution is desired, several steps can 
be taken to accelerate the solution procedure. The details of these steps 
are given in Reference 1. The numerical procedure is based on a point 
-implicit aversion— of the- Jameson, Schmidt and Turkel (Ref 2] finite volume/ 
multistage scheme coupled to a point implicit version of the Jameson 
multigrid method. 

Briefly the idea behind the point implicit method involves modifying 

the original stiff set of equations to another set where the time stiffness 
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associated with the chemical time scales has been removed. This can be 
accomplished by multiplying the unsteady terms by a scaling matrix. For 
example, instead of solving the stiff given by 

3t 3x 3 y + H (1) 

we could solve 

= 3u 3 g 

s T- = - ^7 - + H (2) 

at jX ?Y 

where S is picked to be the- des-ired -scaling matrix. A convenient choice 
for S is the matrix obtained by numerically treating the chemical source 
terms implicitly. Hence the term point implicit. The technique translates 
to solving the equations in pseudo time marching each state quantity at its 
own respective characteristic time scale. Note equations {1} and (2) have 
the same steady state solutions. 

- -Two flame -holder configurations are— considered in this study using 

the compressible' flow equations with H^-Air finite rate chemistry. The 

first uses a wall-generated, oblique shock (ramp) to initiate and hold a 

flame. Examples of this type of reacting flow are shown in Figures 1 and 

2. Figure 1 shows the species contours for a reaction initiated by an 

oblique shock (M. =3.5, HR = 0.026, 4 = 0.1, T. = 90CK, T . = 1Q0CK, with 62'* 

m in lg 

of the available burned) „ Note the flow is entirely supersonic.. If as 
shown in Figure 2 the Mach number is reduced to 2.6 (HR - 0.1, 4=0.1, 89*> 


H" burn) the ~ flow— field takes on a different phenomenology, here a“subsonic 
bubble forms behind the normal Mach stem. The calculations were done inviscidl*. 
as the Reynolds number typical of the ramp flows of interest is lo"'. using 
the point ii-.piicit/rescaling method discussed above, convergence was achieved 
in 200 iterat ons compared to 2 x 10 8 iterations for a purely explicit scheme. 



ff7 - 


3 


The second class of problems involves flame holding behind a rearward 
facing seep. Here the recirculation zone acts to anchor the flame. In this 
case the viscous effects are significant and the full Navier-Stokes equations 
are needed. At the time of the writing of this abstract the rearward facing 
step calculations were not completed, but they would be included in the paper. 
For both flame holder configurations many phenomenologies are expected, and 
an attempt will be made to characterize them. 

Table of Symbols 

HR ’/ Y x HF. / (Cv. X T. + 0.8 x U 2 ) 

L i in m 

l 

<t> Stochiometric mixture ratio 

Y^ P of species i/mixture density 

HF Heat of formation 

■ References 


[I] Bussing, T.R.A . , Murman, E.M., "A finite volume method for the 
calculation of compressible chemically reacting flows," AIAA Paper 85-0331, 
Aerospace Meeting, Reno. 

[Zi Jameson, A., Schmidt, W. , Turkel, E. , "Numerical solution of 
the Euler equations by finite volume methods using Runga-Kutta time-stepping 
schemes," AIAA 81-1259. 

[3] Jameson, A., "Solution of the Euler equations for two-dimensional 
transonic flow by a multigrid method," Princeton University MAE Report No. 
1613, June 198'. 



(if i * i • m; mi (K. <1 turn 


W 


CF v:-- ' 


ricjuro 1 




tin * i i; lin'd ci ii * ci (iiru n<rci cm *cf 






LIBRARY MATERIAL SLIP 

DO NOT REMOVE SLIP FROM MATERIAL 

Delete vour name from this slio when retumina material 

to the library. 

NAME 

DATE 

MS 

Mltlil in IB 


































NASA Langley Form 474 (Rev. Oct. 1999) 











